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Series Solutions of Companding Problems 

By B. F. LOGAN, Jr.* 
(Manuscript received March 10, 1982) 

A formal power series solution (i) x(t) = 2.?m k Xk(t) is given for the 
companding problem (ii) Bf\x(t)\ = my(t), B{x{t)\ = x(t), where B is the 
bandlimiting operator defined by Bg = (Bg)(t) = j- x g(s)[sin X(t - s)]/[ir(t - 
s)]ds and f(t) has a Taylor series with /(0) = 0, /'(0) ^ o. Expressions for the 
x k are given in terms of the coefficients of /, and operations on y, and in a 
different form in terms of the coefficients of the inverse function <f>, <j)\f(x)\ = 
x. A series development is given for a bandlimited z(t), Bz = z, such that the 
solution of (ii) is given by x = B(j>{z). Also a series development is given for 
the "approximate identity", x - B<p{Bf(x)\, where x = x(t), Bx = x, which is 
shown to be a good approximation to x for fairly linear f(x), not necessarily 
having a Taylor series expansion. As an example of one application of the 
results, a few terms are given for correction of the "inband" distortion arising 
in envelope detection of "full-carrier" single -sideband signals. The results 
should prove useful in correcting small distortions in other transmission 
systems. Finally, it is shown that the formal series solution (i) actually 
converges for sufficiently small \m\. This involves proving that the compand- 
ing problem {ii) has a unique solution for arbitrary complex-valued y(t) and 
complex m of sufficiently small magnitude, the solution x{t; m) being, for each 
t, an analytic function of the complex variable m in a neighborhood of the 
origin. It is a curious fact, as shown by an interesting example, that the series 
(i) may converge for values of m for which it is not a solution of (»). 

I. INTRODUCTION 

Suppose x{t) is a bandlimited signal whose Fourier transform van- 
ishes outside the interval [-X, A]. If such a signal is instantaneously 
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distorted by a nonlinear (companding) function f(x), the distorted 
signal f[x(t)\ will, in general, have frequency components outside the 
interval [-X, X]. If the out-of-band components of the distorted signal 
are removed by ideal low-pass filtering, the result is a bandlimited 
signal y(t) whose Fourier transform agrees with that of f\x(t)\ over 
(-X, X). How, and under what conditions, may x(t) be recovered from 
y(t)? When the signals are real-valued, this is known as the compand- 
ing problem of Landau and Miranker (Refs. 1 and 2), hereafter referred 
to as the r.u. companding problem. Before stating their result, and our 
purpose, we introduce some notation. 

The symbol J0 2 (X) will denote the subspace of U. = £*(-», <») whose 
elements are those (square-integrable) functions whose Fourier trans- 
forms vanish outside [-X, X]. Associated with this subspace is the 
bandlimiting operator B x , defined for g in L 2 by 

X, . sin X(t — s) 
. M *(t-«) *■ 

The Fourier transforms of g and Bxg agree over (-X, X), the transform 
of the latter vanishing outside [-X, X]. In the language of Hilbert 
space, Bxg is the projection of g on J^ 2 (X), being the best approximation 
to g in the subspace ^ 2 (X). In case g belongs to ^ 2 (X), we have 

B x g = g. 
It follows that 

Big = B x g, g in U, n = 1, 2, • • • . 

The operator Bx may be applied also to functions belonging to L p , 1 s£ 
p < oo; i.e., to functions g satisfying 

II g \\ P = {J_j£( c ) I Pdt \ < °(1^P< 00 )- 

Here the notation ||g|| p designates the norm of g in L p , or simply the 
Lp-norm of g. The space L x consists of those functions g whose 
magnitude is bounded on the real line, their norm ||g||» being the 
"essential supremum" of \g(t) |, which for functions we will be dealing 
with here, is simply the maximum value of \g(t)\. The operator B x 
may not be applied to an arbitrary bounded function, since the 
associated integral may not converge. However, the integral may 
converge conditionally for a large class of functions; in particular, 
B x g = g, for any constant function g. 
The operator B x is a "contraction" operator on L^; i.e., 

II Bxg II 2^ II* II 2, 
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with equality attaining only for g in ^? 2 (X). This follows from Parsev- 
al's theorem and the definition of B x . 

Applying Schwarz's inequality to the integral equation denoted by 
B\g = g, we obtain the useful inequality 

||£fl.< V\A Wgh, g in m». 
Also, it is easy to show from the integral equation that 
lim g(t) = 0, g in & 2 (\). 

t— »±oc 

We shall also make use of the high-pass operator H x , defined by 

H X = I-B X , 

where J is the identity operator. H x is an identity operator for functions 
h of Z/2 whose Fourier transforms vanish over (—X, X), and it is also a 
contraction operator on L 2 , 

Iltfx£ll2*£ ||g|| 2 , 

with equality attaining only for H x g = g, i.e., for B x g = 0. (In these 
operational equations, is interpreted as the null function.) 

It is clear from the operator definitions and the associated Fourier 
transform relations that any function / in Z^ has the decomposition 

f = g + h, 

where 

g = B x f, h = H x f. 

Since X will be fixed throughout the paper, we will, except where 
emphasis is desired, simply write B, H, and & 2 for B x , H x , & 2 (\), 
respectively. 

Now, using our notation, we may state the important result of 
Landau and Miranker as follows: 

Theorem (Landau and Miranker): Let f(x) be a real-valued function of 
the real variable x, satisfying 
(0 /(0) = 

(ii) < mi s? f'(x) s= m 2 < <», (-00 < x < oo). 
Then to each real- valued y in i^ 2 there corresponds a unique x in & 2 , 
also real-valued, satisfying 

(Hi) Bf(x) = y. 
The solution x of (Hi) may be obtained as the limit of the sequence of 
approximants \x n \ defined iteratively by 

(iv) x n+i = x n - cB\f(x n ) - y\, 
provided only that x t is a real-valued function in & 2 and the real 
constant c is so chosen that 
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(v) max|l-c/'(x)| «r<l. 

X 

The beauty of this result is that, under the hypotheses on /, every 
(r.v.) y in ^ 2 has the representation (Hi) where x is a unique (r.v.) 
function in ^ 2 - In some r.v. companding problems of interest, f(x) 
may not be defined outside some interval and/or the condition on 
f'{x) may not be satisfied over the whole real axis, but rather over 
some interval including the origin. Then the conclusion will apply only 
toy of sufficiently small norm. In such cases, the companding problem 
has two essentially different interpretations. The first is the recovery 
problem: y is known to be of the form (Hi); recover x. The second is 
the design problem: y is a prescribed (desired) signal; find x, if possible, 
so that y is given by (Hi). In this case, one is faced with the problem 
of determining for what y the problem has a solution. 

The speed of covergence of the iterative solution of Landau and 
Miranker is a matter of practical concern. They show that 

l|xn+i - *»ll2 < r||x„ - X n -i\\ 2 . 

Then the constant c in (v) should be chosen to make r as small as 

possible. Assuming that equality may attain on both sides in (»), one 

should choose 

2 . . m 2 -mi 

c = , giving r = ———— • 

mi + m 2 m%-T m\ 

Thus rapid convergence is assured if (m 2 /mi) is not much larger than 
1. If this is not the case, a large number of iterations are, in general, 
required to obtain a close approximation to the solution of the problem. 
In a practical implementation of the iterative scheme of solution (Ref. 
1), the ideal bandlimiting operator is replaced by an approximate 
operator, incurring a certain delay, in addition to (eventually) signifi- 
cant spectral distortions, with the result that the sequence {x„j will 
not converge to the solution x. Thus, in practice, the number of 
iterations to be performed is limited both by practical and theoretical 
considerations. The conclusion is that good approximate solutions to 
companding problems may be conveniently obtained in practice only 
in those cases where the companding function f(x) is fairly linear over 
the range of x(t). 

We should remark at this point that there is only one known 
(nonlinear) r.v. companding problem (see Ref. 3) admitting of an 
explicit noniterative solution; viz., 

B\\og(l + x)} = y, x>-\, x in ^, 

which has a solution if, and only if, the function 1 

+ Here we are applying the B operator to a function not in Z*, the proper interpretation 
being w = 1+ B\-l + exp(-)|. 
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w(t) = B{exp l/2[y(t) + y(t)]\, 

where y is the Hilbert transform of y, extends as a function zero-free 
in the upper half-plane, which will be the case if \\ y\\ 2 is sufficiently 
small. Then the solution is given by 

x(t)=\w(t)\ 2 -l. 

Motivated by the above considerations, pure curiosity, and the fact 
that in many cases of practical interest the companding function and/ 
or its inverse can be well approximated by a polynomial of low degree 
over the range of interest, we are led to consider the case where the 
companding function has a Taylor series expansion, allowing the 
possibility of developing a corresponding series solution to the prob- 
lem. To obtain the terms (1st order, 2nd order, etc.) in the series 
solution it is convenient to multiply y by a scalar parameter m, and 
consider the problem 

Bf(x) = my (1) 

to be solved for x in 3k, given y in ^, for companding functions 

f(x) = lb k x h , \x\<R (2) 

i 

6i * 0. 

For sufficiently small \x\, f will have an inverse <£, 

x = <P\f(x)\ 

<t>(y) = I a k y k , lyl^-Ro*. (3) 

i 

We assume that the solution x = x(t; m) of (1) has a series expansion 
in the parameter m, 

x(t; m) = £ m k x k (t), (4) 

i 

where the Xk(t), aptly described as fcth order corrections, (not to be 
confused with the Landau-Miranker approximants) depend only on 
y(t) and /. Presumably, in cases of small distortion, a few terms of the 
series would give a satisfactory approximation to the solution. 

Explicit expressions for the first five of the Xk(t) are given in the 
sequel, first in formulas involving the coefficients of <f>, and next, the 
coefficients of /, together with certain operations on y. These formulas 
reveal how the Fourier transforms of the Xk(t) may be calculated from 
the Fourier transform of y(t), if this be given. 

Next, we find a series development of z in ^, 
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z(t) = z(t; m) - 2 m*2 fc (t) (5) 

i 

such that the solution to (1) is given (presumably for sufficiently small 
\m\)by 

x = B4(z). (6) 

We find that z x =y,z 2 = 0, and in case a 2 in (3), or b 2 in (2), vanishes, 
we have, in addition, z 3 = z 4 = 0. That is, under certain conditions, 
z = my, implying that B<f>\Bf(x)} is an "approximate identity" (=x) for 
x in ^ 2 , especially if / is odd and fairly linear, or if x{t) is a predomi- 
nantly low-frequency function. 

To further investigate the approximate identity, we introduce the 
parameter m again, and obtain expressions for Uk in 

B4>{Bf(mx)\ = £ mV x in 9k- (?) 

i 

To see how interchanging / and <t> affects the approximate identity, we 
compare u* with i>* in 

Bf\B<t>(mx)\ = £ m k v k , x in 3 % . (8) 

i 

As expected from the series development of z, we find u x = Vi = x, and 
u 2 = V2 = 0. Further comparisons [with the same m in (7) and (8)] 
should be made for the case /'(0) = 0'(O) = 1. For this case, we find 
U3 = U3 = 2b\B(x-Hx 2 ), which may be small if b 2 is small or if Hx 2 is 
small. In case b 2 = 0, we find u k = v k = for k = 2, 3, 4, and u 6 = 
v 5 = 3blB(x 2 .Hx 3 ). 

These series developments of the approximate identity suggest that 
it would be useful in obtaining an approximate solution to the r.v. 
companding problem for fairly linear companding functions, not 
necessarily having a Taylor series expansion, but merely satisfying 
/(0) = and 

< m 1 ss f'(x) m 2 < oo, (-00 < x < oo). (9) 

Compelled by this suggestion, we digress in the Appendix to show for 
such / that 

\\x-B<t>lBf(x)\\\ 2 <y\\x\\ 2 , x in J%, (10) 



where 



e 2 m 2 



7 4(1 + «)' 6 rn r 



(Note that y = 1/8 for m 2 lm x = 2.) Thus in many companding 
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problems, B(f>(y), involving only one filtering operation, would be an 
adequate approximation to the solution x. We go on to define an 
iterative procedure, involving both / and its inverse 0, obtaining 
approximants converging to x for y < 1, offering an alternative to the 
solution of Landau and Miranker in cases where (m 2 /mi) < 3 + 2\/2. 
In any case, B(f>(y) is suggested as a good choice for X\ in their iterative 
solution. We note, in leaving this topic, that the inequality (10) is 
invariant to the interchange of / and its inverse <f>. 

Returning to the series solution, we apply the results to the problem 
of compatible single-sideband transmission (Ref. 4), obtaining a few 
terms for correction of the "in-band" distortion arising in envelope 
detection of "full-carrier" single -sideband signals. 

Although the original intent of the work here was to obtain expres- 
sions for the first few terms of the series solution (4), supposedly 
adequate for correcting small distortions, the mathematical question 
naturally arises in the end as to whether the series actually converges 
for sufficiently small \m\ (or equivalently, for \m\ = 1 and ||y|| 2 
sufficiently small), or whether it is merely an asymptotic series. It is 
indeed a pertinent mathematical question, since the expressions for 
the Xk(t) were obtained by purely formal manipulations of power series 
and application of the operators B and H. The resulting expressions 
become progressively cumbersome and complicated, with no obvious 
general form, offering no possibility of establishing (from them) 
bounds on | Xk(t) | which would ensure convergence of the series. The 
remainder of the paper is addressed to the problem of establishing the 
convergence of the series. 

If we suppose in the r.v. companding problem that the series con- 
verges for real-valued m of sufficiently small magnitude, then it would 
also converge for similar complex m, suggesting that the companding 
problem [for fixed y(t)] would have a solution for all complex m of 
sufficiently small magnitude. This, in turn, suggests that the problem 
would have a solution for arbitrary complex-valued y(t) in ^ 2 and all 
complex m of sufficiently small magnitude, depending on / and the 
norm of y. That this is a fact has been established previously (Ref. 3) 
only for complex- valued y(t) whose Fourier transforms vanish outside 
[0, X], (or [—X, 0]), the Fourier transform of the solution x(t) having 
the same property. In this case (with m = 1), the solution is given by 
x = B<t>{y) for y of sufficiently small norm; i.e., in case the Fourier 
transform of x(t) vanishes outside [0, X], (or [—X, 0]) the "approximate 
identity" is an exact identity, 

x = B<f>[Bf(x)} = <£{/(*)} 

for x of sufficiently small norm. This result can be explained, roughly, 
by the fact that nonlinear (analytic) distortion of such x{t) does not 
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produce both "sum and difference" frequency components, but only 
"sum" components. 

In order to prove that the series solution actually converges for 
sufficiently small |m|, we show that the companding problem (1), 
where y(t) is an arbitrary complex-valued function in ^ 2 , has a solution 
x(t; m) for all complex m of sufficiently small magnitude, this solution 
being, for each fixed t, an anaytic function of the complex variable m, 
from which it follows that the solution has a Taylor series expansion 
in m; i.e., 

x(t;m) = J m k x k (t), \m\ < m. (11) 

i 

To obtain this result, we first have to establish for the complex- 
valued (c.v.) companding problem the analogue of A. Beurling's 
uniqueness theorem (see Ref. 1) for the r.v. companding problem. 

We then examine in detail a specific problem illustrative of the 
theory and some of its nuances; viz., the problem (taking X = 2 for 
convenience) 

« I X I sin It /i0 v 

B W - m ~2T <12) 

for which the solution is (at least for sufficiently small \m\) 

nn sin 2t ., / sm A t*n\ 

x = x(t;m) = 2/3 -^- - /3 2 I — I , (13) 

where 

/? = m/(2 + m). 

The rather surprising revelation of this example is that, although the 
series expansion in m of x(t;m) converges, uniformly in t, for | m \ < 
2, it is not a solution of (12) for all such m. Furthermore, one might 
reasonably assume that (13) is a solution of (12) for all m other than 
-2, but this is not true either. As an illuminating exercise, we deter- 
mine precisely the set of m for which (13) is a solution of (12). 

II. THE INVERSE SERIES METHOD 

To obtain a series solution to (1), we first think of recovering from 
y(t) the out-of-band components of f\x(t)\, so that we might apply the 
inverse function <f> to the whole in order to recover x(t) = x(t;m). We 
have 

f\x(t)) = my(t) + h(t), (14) 
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where h(t) = h(t;m) is some unknown "high-pass" function satisfying 

Bh(t) = (15) 

and hence 

x(t) = <t>\my(t) + h{t)\. (16) 

It is convenient at this point to introduce the high-pass operator 
defined by 

H = I-B, (17) 

where J is the identity operator. Thus applying H to (16) we have 

H<t>\my(t) + h(t)\ = 0. (18) 

We would like to solve (18) for h(t), which we think of as small in 
cases of interest. 
Now we assume that 

0(y) = £ a*y* for sufficiently small \y\, (19) 

i 

a x * 
and that in (16) 

h(t) = h(t;m) = £ m k h k (t), Hh k = h k , k& 2, (20) 

2 

x(t) = x(t;m) = £ m k x k (t), Bx k = x k , k& 1, (21) 

i 

where hk(t) and x k (t) do not depend on m. 

We want to expand <f>\my{t) + h{t)\ as a power series in m. To do 
this it is convenient to write 

my(t) + h(t) = £ m k h k (t), (22) 

i 

where we identify 

y(t) = h x (t) in @2. (23) 

Then we write 

\ S m k h k (t) \ = F(m;t) = £ m k F k (t). (24) 

For convenience we suppress the variable t and write simply x k , h k , F k . 
In terms of the coefficients a* in 

COMPANDING PROBLEMS 3015 



0(y) = 2 w k , 

i 

we find, equating coefficients of m k in (24), 

JPi = a^ (25.1) 

F 2 = ai/ia + <M (25.2) 

i? 3 = ai h 3 + a 2 (2h 1 h 2 ) + cM (25-3) 

F 4 = a x h 4 + a 2 (2hji 3 + hi) + a 3 (3/if^) + a 4 h\ (25.4) 

*6 = diha + a 2 (2hih 4 + 2h 2 h 3 ) + a 3 (Sh 2 h 3 + 3/ii/ii) 

+ a 4 (4fc?/i 2 ) + a 5 fc? (25.5) 

F 6 = a x /i 6 + a 2 (2/ii/i 5 + 2/i 2 fc 4 + ^1) + a 3 (3/i?ft 4 + ShJhfa 

+ hi) + a 4 (4h\h 3 + 6/if/ii) + ab(6MAi) + aM. (25.6) 

F 7 = ai/i 7 + a 2 (2/ii/i 6 + 2fc 2 h 5 + 2ft 3 h 4 ) 

+ a 3 (3/i?/i 5 + 6/ii^/i 4 + 3/iifti + 3fti/i 3 ) 

+ a 4 {4h\h 4 + 12hlh 2 h3 + 4hM) 

+ a 5 (5h\h 3 + 10hlh 2 2 ) + ae(6hlh 2 ) + a n h\ (25.7) 

F s = aih s + a 2 (2/iih 7 + 2/i 2 /ie + 2h 3 h 5 + hi) 

+ a 3 (3fcf/i 6 + 6/ii/i 2 /i 5 + 6h 1 h 3 h 4 + Shlh 4 + Sh 2 hl) 

+ a 4 (4h\h 5 + I2h\h 2 h 4 + 6hlhl + I2hj%h 3 + K) 

+ a 5 (5h$h 4 + 20hih 2 h 3 + 10h 2 M) 

+ a 6 (6hlh 3 + 15/if/ii) + (nOhfh,) + a B hl (25.8) 

F 9 = aih*i + a 2 (2hih 8 + 2h 2 h n + 2h 3 h 6 + 2h 4 h 5 ) 

+ a 3 {Shlh 7 + eh^he + 6hth 3 h 5 + Zh x h 2 4 + 3hlh 

+ eh 2 h 3 h 4 + h 3 3 ) + a 4 (4hih 6 + 12h?/Wi 5 + 12h\h 3 h 4 

+ 12hMh 4 + 12/i a /2 2 /i§ + 4hlh 3 ) + a b (bh\h 5 

+ 20h\h 2 h 4 + lOhlhl + 2,0h\hlh 3 + bh^) + ae(6hlh 4 

+ 30hih 2 h 3 + 20h%hl) + a 7 (7/if/i 3 + 21h\hl) 

+ a 8 (8hlh 2 ) + a*hi (25.9) 

Fw = a x hiQ + a 2 (2/ii/i9 + 2h 2 h& + 2h 3 h 7 + 2h 4 h 6 + hi) 

+ a 3 (Shlh 8 + §hih 2 h n + 6hih 3 he + 6h x h 4 h 5 + 3h!/i 6 
+ 6h 2 h 3 h + 3/i 2 /i 4 + Sh 2 3 h 4 ) + a 4 (4/i?^7 + 12hlh 2 he 
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+ 12/i?/i 3 /i 6 + 6h 2 h 2 4 + \2hMh 5 + 2Ah x h 2 h 3 h t + 4h t hl 

+ 4h 3 2 h 4 + 6h 2 2 h%) + th(5hthe + 20h\h 2 h b + 20h\h 3 h< 

+ Z0h 2 Mh A + 30h 2 h 2 hl + 2QhMh 3 + hi) + a 6 (6/i?/i 5 

+ SOhihifa + IShihl + GQhlhZka + 15h 2 h$) 

+ ch{lh\h A + A2h\h 2 h z + 35/#ii) + a 8 (8h\h 3 

+ 28/i?/i|) + 09(9/1?^) + a 10 h\°. (25.10) 

In general the coefficient of a m in the expansion of F n consists of sums 
of products of the hk corresponding to partitions of n into m parts. 
The coefficient of the product is m! divided by the product of the 
factorials of the exponents of the hk (the multinomial theorem). For 
example, in Fi the coefficient of a 5 is found by writing down the 
partitions of 10 into 5 parts and proceeding thus (see Table 24.2, Ref. 
5): 

l 4 , 6 -> hit* coef. = 5!/4! = 5 



l 3 , 2, 5 


hihdn*, 


5!/3! = 20 


l 3 , 3, 4 


hifi3h,4 


5!/3! = 20 


I 2 , 2 2 , 4 


hih 2 h A 


5!/2!2! = 30 


I 2 , 2, 3 2 


hih 2 h3 


5!/2!2! = 30 


1, 2 3 , 3 


hih 2 h 3 


5!/3! = 20 


2 6 


h\ 


5!/5! = 1. 



Now we may obtain a formal series solution (21) by successively 
solving for the hk by requiring 

HF k (t) = 0, k = 1, 2, • • • (26) 

and setting 

x k (t) = BF k (t) = F k (t), k = 1, 2, • • • . (27) 

Recall that hi = y, the given bandlimited (low-pass) function, and all 
the other h k are high-pass functions. 1 We have Hh k = h k , k 2* 2 and 

HFi = diHh! = (28.1) 

HF 2 = ai h 2 + a 2 Hh\ = 

h 2 = - ^ Hh\ (28.2) 



* Actually, for k ^ 2, hi, is a bandpass function whose Fourier transform vanishes 
over (-X, X) and outside [-kX, kX]. This can be seen from (28.1)-(28.k). 
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HF 3 = aJh + 2azH(h 1 h 2 ) + aM. = 

h 3 = H Hihx.Hti) - ~ Hh\ (28.3) 

a{ ax 

HF 4 = aihi + 02{ZH{h x hz) + Hhl) + daaHihlht) + a<Hh\ = 
h 4 = -^4 H[hi-H(hi.Htib] + ^r BUh'Mb 

- feY H[Wi) 2 ] + ^ht H(«-flM) - - ^M (28-4) 
\ai/ a x a x 

HF 6 = d/iB + (h[2H(hM + 2H(h 2 h 3 )] + a 3 [3H(hlh 3 ) 

+ ZH{h x h%)] + 4a A H{hlhi) + aM = 

at ai 

+ ^ H[hr-Wl) 2 ] - *§* mhi-Hiti-Hhl)] 

CLx All 

+ ^ fl(At.flM) + ^r mHh\)-H{h x .Hh\)] 
ai a x 

- ^ H[W 2 )-Hhl] - ^ HlhlHih^Hhl)} 

ai fli 

ai ai 

+ ifW4 H{h i. Hh i) _ 55 flfcf. (28.5) 

af a x 

Now replacing /ii by y we have from (25), (27), and (28) 

X, = BF l = a x y (29.1) 

x 2 = BF 2 = a 2 By 2 (29.2) 

x 3 = BF 3 = 2a 2 B{yh 2 ) + a 3 By 3 

x 4 = BF 4 = a 2 (2B(yh 3 ) + Bh 2 2 ) + Za*B(y 2 hi) + a 4 By 4 

= t4B[y.my.Hy*)}-^B(y.Hy>) 
af a x 
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+ jj B(Hyr - ^ B(*W) + a t By> (29.4) 

x b = BF 5 = a 2 [2B(yh 4 ) + 2B(/i 2 /i 3 )] + 3a[B(y 2 h 3 ) + B(yh 2 2 )] 
+ 4a<B{y*h 2 ) + a^By 6 

= -^i B{y.H[y-H(y-Hy 2 )]\ + ^ B[y.H(y.Hy 3 )] 
- ^ B[y.H(Hy*?] + ^ B[y . H( y jy,, 

"1 Ol 

+ ^ Bl(Hy>).(Hy>)) + «2b B[y^H(y.Hy')) 

_22l s^.^3) + **S B[y.(Hy^] 
«i ai 

4a 2 a 



d 



-B(y 3 .//y 2 ) + a5By 5 . (29.5) 



If in (29) we replace the H operator by / — B and collect terms we 
obtain 

xi = a x y (30.1) 

x 2 = azBy 2 (30.2) 

x 3 = ^ B^By 2 ) + fa 3 - ^) By 3 (30.3) 

{ 2 nr _. „,.. ,,.2x1 / 4fl i 2a 2 a 3 



x 4 = -r B^-Bdy-By 2 )] - -f - == B(y.By 3 ) 
d\ \ai ai 

+ | B(By°) 2 - M - ***) B(y*.By*) 
a x \ai ai / 

(5a 2 a 3 5a 2 \ ^ . 
a 4 — + ~r By 4 30.4) 

a, af/ 

« = ^ B | y .fli y .B(y.«,«)]| + fe - *4) Bly.Biy.By>)] 
/6^_1^\ [y ^ W)] + 2^ 



\ af a? / ~ l " ~" " J n ' a? 
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\ o,\ a x ] 

+ /2fiB. _ lOdfi + ««\ B( B/) + 1^ B[(Bf) 

\ a l a? ai / ai 



/6ala3_12a2\ „, , 0/ „ ^ 
V a? a\ 



■B(y-By 2 )] + ^ - ^ p3[y 2 -B(y-By 2 )] 



+ 



/2a|j _ 4^\ (By3)] 

\fli ai/ 



+ /3a| _ 12^ + 12a| \ 

V ai ai ai J 

+ /4^_20a|a3 + 20a|\ B(y 3. B3/2) 
\ a x a\ af J 

+ /«sa + S*_M-!_a4\ By6 . (30 . 5) 

\ ai af a? ai / 

Note that if y belongs to MX/n), then By k = y h for k = 1, 2, • • • , n. 
In this case we will have x n = a„y n . So the sum of all coefficients in 
the expressions for x n must be a n . If </> is an odd function these formulas 
simplify considerably. It is rather curious that if a 2 = 0, a 3 # 0, the 
coefficient of B{y 3 -By 2 ) vanishes, whereas the coefficient of B(y 2 - 
By 3 ) does not. The coefficients in (30) are more simply expressed in 
terms of the coefficients in the power series for / as we see below. 

III. FORWARD SERIES METHOD 

We can also solve (1) in the "forward" direction by writing 

Bf(mxi + m 2 x 2 + m 3 x 3 + • • • ) = my, (31) 

where 

fix) = l b& k . (32) 

i 

Then applying the expansion (24) and (25) to /(£ m k x k ) we have, 
equating coefficients of m k , 

Bb lXl = y (33.1) 

Bb lX2 + Bb 2 x\ = (33.2) 

Bb,x 3 + Bb 2 (2 Xl x 2 ) + Bbaxl = (33.3) 

Bb lX4 + B&2(2xiX3 + x\) + Bb^x\x 2 ) + Bb<x\ = (33.4) 

3020 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1983 



Bbix s = Bb 2 (2 Xl x 4 + 2x 2 x 3 ) + Bb 3 (3xlx 3 + 3x lX 2 2 ) 

+ Bb 4 (Ax\x 2 ) + Bb b x\ = 0. (33.5) 
Then solving (33) successively for x k , (Bx k = x k ), we have 

xi = y/bi (34.1) 

*2 - - p B y 2 (34.2) 

x^^Biy-By^-^By* (34.3) 



h *B[yB(yBy*)) + ^* - • 



x 4 = --f B[y.B(y.By 2 )) + -^ B(y-By*) 



-%B { Byr + SJ ^B(y>.By>)- b -±By< 

0\ Oi Oi 

x 5 = ^B\y.B[y.B(yBy 2 ))\ + ^ B[(By 2 ).B(y.By 2 )] 
+ ^ B[y.B«yV| - ** S[yB(y5y 3 )] 
- S^ B[y.5(y 2 .fiy 2 )] - ^ B[(By 2 ) . (By 8 )] 



B[y 2 -5(y-fly 2 )]--JrB[y.(fly 2 ) : 



(34.4) 



6! l " w " /J 6! 
+ ^B(y.By )+ ^B(y3. By2) 

+ ^B(y 2 -By 3 )-^By 5 . (34.5) 

Oi Oi 

These correspond to the solutions for the hk in (28) with H and B, and 
a's and 6's interchanged, except here X\ = y/6 x as compared by h\ = y 
in (28). They agree with the formulas in (30) according to the identities 
in reversion of series (see 3.6.25, Ref. 5) 

flioi = 1 (35.1) 

a\b 2 = -a 2 (35.2) 

a% = 2a| - aia 3 (35.3) 

a\b 4 = 5a 1 a 2 a.' t - a\a 4 — 5a 2 (35.4) 
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a?o 5 = 6a?a 2 a 4 + 3afa§ + Ua 2 4 - afas - 21a 1 aia 3 (35.5) 

a"68 = 7aia205 + 7a 3 a 3 a 4 + 84aia 3 a 3 - a\ae 

- 28a?a 2 a§ - 42ai - 28a?a|a 4 (35.6) 
a} 3 6 7 = 8ala 2 a« + 8a?a 3 a 5 + 120a?aia 4 

+ 180a?aiai + 132a| - afa 7 - 36a?a|a 5 

- 72a?a 2 a 3 a 4 - 12a?a 2 a 3 a 4 - I2a\al - 330aiafa 3 . (35.7) 

Of course the a's and 6's may be interchanged in (35). Actually the 
expressions for x k in (34) are analogous in a way to the coefficients a k 
expressed in terms of the b k as given by (35) (with the a's and 6's 
interchanged). That is, if it were not for the B operator in eq. (33.k), 
we would have simply x k = a k y k according to the determining equations 
for the inverse coefficients. Because of the B operator, successive 
solutions for the x k generate powers of y interposed with B operators 
and powers of \B( ■ )\ in all combinations so that, for example, in x 5 we 
have a number of terms with coefficients o 2 6 3 /6? that combine only 
when B is replaced by J to give -216 2 6 3 y 5 /6i» corresponding to the last 
term in (35.5). Note if all the b k = 1, the sum of the integer coefficients 
in x n is (-l) n+1 as this is the case y = fix) = x/(l - x), x = 0(y) = y/ 
(1 + y). Note also that in the expression (34.5) for x 6 , for example, 
there are five groups of functions with common "6" coefficients. These 
combine with certain weights, depending only on the coefficients of 
f(x), to give x b . Similarly, in the expression (29.5) obtained from the 
inverse function, there are again five groups of functions with common 
"a" coefficients that combine with certain weights, still depending 
only on the coefficients of f{x), to give x b . The interesting and rather 
puzzling fact is that the groups are not identical but overlap. 

IV. A SOLUTION OF THE FORM x(t) = B<t>[z(t)) f z in J0 2 

For the solution to (1) we have x = £?° m h x k , where according to 
(25) and (29), we have 

Xi = a x y 

x 2 = ai/12 + a*y 2 = (hBy 2 

x 3 = ai/i 3 + 2a 2 y/i2 + a 3 y 3 = laTPiyh*) + a^By 3 . 

Since 

0(my) = ai/ny + a^m 2 y 2 + a 3 m 3 y 3 + • • • , 

we see that x = B<t>(my) + J0{m 3 ), m -+ 0. Then setting m = 1 (with y 
sufficiently small) we can conclude that 
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x = B<f>(y) + My\ y->0. (36) 

At least as y — > 0, B$(y) is a better approximation than <f>{y) = x + 
J?(y 2 ). Also B<j>(y) could be a good approximation to x without y being 
small, as would be the case if y were a predominately low-frequency 
function (compared with its top frequency). This suggests that given 
y in (1) we determine a bandlimited function z, perhaps close to y, 
such that the solution to (1) is given by 

x = B<f>(z). (37) 

To determine a series solution for z we set 

00 

z = £ z k m\ Bz k = z k (38) 

i 

and expand </>(z) in a power series in m: 

timz! + m 2 z 2 + m 3 z 3 + •••) = ! m h F k . (39) 

1 

The Fk are given by (25) with z* replacing /i*. The difference now is 
that the F k are bandlimited to [— k\, kX], i.e., <t>(z) is not bandlimited 
(in general) with z in i^ 2 . However, we must have 

BF k = x ht (40) 

where in terms of operations on y the Xh are given conveniently by 
(29). We have 

flctiZi = xi = axy (41.1) 

B{a x z 2 + a 2 z 2 ) = x 2 = a^y 2 (41.2) 

o fl 2 

B(a lZs + 2a 2 z l22 + a 3 zl) = x 3 = B(y-Hy 2 ) + a^By 3 (41.3) 

B[aiz 4 + a 2 (2ziz 3 + zl) + a 3 (3z?z 2 ) + a 4 zi] = x A 
.*§BLy.H(y.Hy*)\- — B{y.Hy*) 

+ a \ B(Hy 2 ) 2 - 5^3 B{y K H y*) + a 4 fly 4 (41.4) 

af ai 

B[a x z b + a 2 (2ziz 4 + 2z 2 z 3 ) + a 3 (3ziz 3 + 3z x z 2 ) 

8a 4 
+ a 4 (4z?z 2 ) + a^] = x 5 = - -£ B{y.H[y.H(y.Hy 2 )]) 

+ Up Bly-myHy*)) - ^ B[ y .H(HyV] 
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+ ^ B[yH(y*.Hy*)] - ^ B(yHy') 
ai a x 

- ^ B[(Hy 2 ).H(y.Hy 2 )) + ^ B[(Hy 2 ) ■ (Hy*)] 

+ ^ B[y*-H{yHy*)] - ^ B(y*-Hy°) 

+ ^ B[y.(Hy 2 ) 2 ] - ^ B(y 3 .Hy 2 ) + a 5 By 5 . (41.5) 

Solving these equations successively for z* we find 

Zl =y (42.1) 

2 2 = (42.2) 

z 3 = -^iB(y.Hy 2 ) (42.3) 



* = *4 B[yH(yHy>)) - ^ B(yHy 3 ) + $ B(Hy > 



2\2 



- 2^ B{y*-Hy*) + ^ B[y.B(y.flSy 2 )]. (42.4) 

4ai 
The first and last terms in (42.4) combine (H + B = I) to give — 

fli 

B(;y 2 -.H3' 2 ), which then combines with the fourth term to give 

24 = ^| (4ai - Sa iaa )B(y 2 -Hy 2 ) - ^ B(y.Hy>) 
a.\ ai 

+ ^| B(Hy 2 ) 2 (42.4a) 

fli 

e 5 = -^fi{yH[3'-H(y.Hy 2 )]} 
at 

- ^ B[y.B(y 2 -i/y 2 )] + - 3 (4a 2 - - Sa^B^-Hy 3 ) 
at ai 

+ ^(3 ai a 3 -2a 2 2 )B[y.(Hy 2 ) 2 ) 
at 

+ ^ (3a 2 a 3 - ai a 4 )fi(y 3 .i/y 2 ) - ^ B(y.Hy*) 
a\ fli 
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- ^4 B[(Hy 2 ).H(y.Hy 2 )] + ^ B[(Hy 2 ) . (Hy 3 )]. (42.5) 

The expressions for the third- and fourth-order terms, z 3 and z 4 , are 
quite simple, owing to the fact that Z\ = y, Z2 = 0. Note in eq. (41.n) 
the "diagonal" term a„Bz" on the left is cancelled by the term a n By n 
appearing in x n on the right. Also in case o 2 = we have 

2i - y (43.1) 

z 2 = (43.2) 

23 = (43.3) 

z 4 = (43.4) 

z 5 = -^B(y 2 .Hy 3 ). (43.5) 

V. THE APPROXIMATE IDENTITY 

The series development in the previous section suggests that as a 
practical expedient one might take 

x = B<f>(y) = BMBf(x)} = <p[f(x)\ = x. 

That is, what appears to be the naive thing to do may in fact be quite 
good, especially for odd functions <p (or /) that are not severely 
nonlinear. The interposition of the bandlimiting operator between a 
nonlinear function and its inverse and then subsequent bandlimiting 
is an interesting "approximate identity" that we examine further in 
the Appendix. One might ask how the interchange in the order of a 
particular function and its inverse in the transformation affects the 
approximate identity. The series expansion of the approximate iden- 
tity may shed some light on the general problem. To keep track of the 
various orders, it is convenient to introduce the parameter m as before. 
We have 

Bf(mx) = J m k bkBx k , (44) 

i 

<t>\Bf(mx)\ = i m"F k , (45) 

i 

where the F k are given by (25) with b k Bx k replacing h k , and Bx = x. 
Thus, 

F l = aA* = x (46.1) 
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F 2 = ai b 2 Bx 2 + a 2 b\x 2 (46.2) 

F 3 = atbaBx* + 2a 2 b 1 b 2 x-Bx 2 + a 3 b\x 3 (46.3) 

F 4 = Oi&iBx 4 + 2a 2 b 1 b 3 x-Bx 3 + thb&Bx 2 )* 

+ Za 3 b 2 b 2 x 2 .Bx 2 + a A b\x A (46.4) 

F 5 = a l b b Bx b + 2a 2 b x b 4 x-Bx 4 + 2a 2 b 2 b 3 (Bx 2 ) ■ (Bx 3 ) 

+ Za 3 b\b 3 x 2 -Bx 3 + Sa^bjx-iBx 2 ) 2 

+ Aa 4 b\b 2 x 3 -Bx 2 + a b b\x\ (46.5) 



Now we set 



where 



B<f>\Bf(mx)\ = £ m k u k , (47) 

i 



u k = BF k . (48) 

Now note in (46) that if all the terms in BF k involving x were of the 
form Bx k then BF k would vanish identically for k 2s 2 because tf>{/(x)} 
= x. So we will introduce the high-pass operator H = I - B to collect 
the terms Bx k that cancel. For example, to collect terms fix 3 in BF 3 
we write 

B{x-Bx 2 ) = B(x 3 - x-Hx 2 ) = Bx 3 - B(x-Hx 2 ). 

Thus, 

Ul =BFl = x (49.1) 

U2 = BF 2 = a x b 2 Bx 2 + a 2 b\Bx 2 = Q (49.2) 

W3 = sf 3 = aiboBx 3 + 2a 2 b 1 b 2 B(x-Bx 2 ) + a 3 b\Bx 3 
= axbsBx 3 + 2a 2 b 1 b 2 B(x 3 - x-Hx 2 ) + a 3 b\Bx 3 
= -2a 2 b x b 2 B{x-Hx 2 ) (49.3) 

u 4 = BF 4 = a x b 4 Bx 4 + 2a 2 b 1 b 3 B(x-Bx 3 ) + a 2 b\B{Bx 2 ) 2 
+ 3a 3 b 2 b 2 B(x 2 -Bx 2 ) + a 4 b\Bx 4 
= ai b 4 Bx 4 + 2a 2 b 1 b 3 B(x 4 - x-Hx 3 ) + a 2 b\B(x 2 - Hx 2 ) 2 

+ Sa 3 blb 2 B(x 4 - x 2 -Hx 2 ) + a 4 b\Bx 4 
= -2a 2 b x b 3 B(x-Hx 3 ) - (2a 2 b 2 2 + 3a 3 b 2 b 2 )B(x 2 ■ Hx 2 ) 

+ a^lBiHx 2 ) 2 (49.4) 
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u B = BF 5 = ai b 5 Bx 5 + 2a 2 b 1 b 4 B(x-Bx 4 ) 
+ 2a 2 b 2 b 3 B[(Bx 2 )-(Bx 3 )] 
+ 3a 3 b 2 b 3 B(x 2 -Bx 3 ) + 3a 3 6i&!B[x.(£x 2 ) 2 ] 
+ Aa 4 b\b 2 B(x 3 -Bx 2 ) + a b b\Bx b 
= a x b b Bx b + 2a 2 b 1 b 4 B(x 5 - x-Hx 4 ) 

+ 2a 2 b 2 b 3 B[(x 2 - Hx 2 ).(x 3 - Hx 3 )] 
+ 3a 3 b\b 3 B(x b - x 2 -Hx 3 ) 

+ Za 3 bMB[x-{x 2 - Hx 2 ) 2 ] + 4a 4 blb 2 B(x 6 - x 3 -Hx 2 ) 
+ a b b\Bx b 
= -2a 2 6 1 fe 4 5(x.i/jc 4 ) - (2a 2 b 2 b 3 + %a 3 b\bz)B{x 2 .Hx 3 ) 
- (2a 2 b 2 b 3 + 6a 3 6 x 6l + Aa 4 b\b 2 )B(x 3 -Hx 2 ) 
+ 2a 2 b 2 b 3 B[(Hx 2 ) . (Hx 3 )] + Za 3 bMB[x- (Hx 2 ) 2 ]. (49.5) 

Now in order to assess the symmetry or lack of symmetry in inter- 
changing <j> and / we can use the identities (35) to express the mixed 
coefficients of u n in terms of the b k or the a k alone. We have 



(50.1) 
(50.2) 



Ui = x 

1/2 = 

u 3 = 3 dB(x-Hx 2 ) (50 3 ) 

u 4 = 4 C 1 B(x-Hx 3 ) + 4 C 2 B(x 2 -Hx 2 ) + .CMHx 2 ) 2 (50.4) 
u 5 = sC^x-Hx 4 ) + b C 2 B(x 2 -Hx 3 ) + b C 3 B(x 3 -Hx 2 ) 

+ 5 C 4 B[(Hx 2 ) ■ (Hx 3 )] + 5 C 5 B[x. (Hx 2 ) 2 ], (50.5) 

where 

3 Cj - -2026x62 _ — = — (50 3a) 

n _ o k k 2 ^ 2 ^3 2a 2 a 3 4a| 
4 C X - -2036x63 = -75- = —t- - -x (50.4a) 

U\ U.\ U\ 

4 C 2 = -(2a 2 b\ 4- 3a 3 b 2 b 2 ) = -^ + ^l = _M + ^l 



bl "6? at T a? 



.3 



bl a? 



n — « 1.2 _ w 2 02 

4^ 3 — 0202 7a = ~ 
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c -_ 2aD0 _2^ = _10o|o3 2a^4 10a| ( } 

46|b 3 361 

5 C 2 = -(2a 2 6 2 6 3 + 3a 3 6?6 3 ) = - —gr + ~yT 

= (2al - aia 3 )(2ai - 3a!a 3 )/a? 

5 C 3 = -(2a 2 6 3 + 6a 3 6ife 2 + 4a 4 6f6 2 ) 

_ 86| _ 126j6 3 46 2 6 4 = 4aj _ 8aia 3 4a 2 a 4 
~"&f 6? + 6f ol ol a? 

26j6 3 _ 2aja 3 _ 4a| 
5C4 = 2a 2 6 2 6 3 — — , 3 — 7 « 

Oi cZi "1 

6 °2 3616 3 _ 3a|a 3 
5C5 — 3a 3 6i6 2 — , 4 . 3 — 7 ■ 

Ol Oi ai 

Now if we interchange the order of <t> and / in (45) and write 

Bf\B4>(mx)} = X m k Vk, ( 51 ) 

1 

we obtain the u k by replacing u k by v k in (49.k) and then interchanging 
the a's and 6's. We have u x - «i - x and u 2 = o 2 = 0. We should 
compare u k and v k , k&3, for /'(0) = 1 = d = 6 X . Then we have 

Uft . o, ft = 1, 2, 3. (52) 

But we have, for example (with ai = 61 = 1), 
u 4 = 2b 2 bsB(x.Hx 3 ) + (36 2 6 3 - 46 3 2 )B(x 2 • Hx 2 ) - b\B{Hx 2 ) 2 (53.1) 
v 4 = (26 2 6 3 - 4bl)B(x-Hx 3 ) 

+ (36 2 6 3 - 2bl)B(x 2 .Hx 2 ) + b\B(Hx 2 ) 2 . (53.2) 

If, however, 6 2 = (a 2 = 0) we have 

u k =v k = Q, k = 2, 3, 4, (53.3) 

and if d = 61 = 1, 

u 5 = Vs = SblB(x 2 .Hx s ). (53.4) 

In case 6 2 = 64 = 6 6 = 0, (a 2 = a 4 = a 6 = 0), we have 

1h 3 

-*gB[x.{Hx*n (54) 

0l 
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where the coefficients expressed in terms of the a's are 





36365 
b\ 


3a 3 a 5 
a! 


9a! 
a?' 






56365 
6? 


96| 

6? " 


50305 
a? 


6a! 
a? 


36! 
6? " 


3a! 
" a? 



So we do not have, in general, u 7 = v-j for odd functions /(x) with 
/'(0) = 1. 

VI. APPLICATION TO COMPATIBLE SINGLE-SIDEBAND TRANSMISSION 

The mathematical problem of compatible single-sideband transmis- 
sion was formulated in Ref. 4. Given a signal y(t) in ^ the problem 
is to determine m such that the equation 



B{V(1 + s{t)) 2 + s(t) 2 \ = my(t) + 1 (55) 

has a solution s(t) in ^ 2 - In (55) s(t), sometimes called the quadrature 
signal, is the Hilbert transform of s(t), and B is operating on the 
envelope of the single-sideband signal. The idea is to transmit a single- 
sideband signal that is compatible with receivers designed for double- 
sideband (AM) reception. Setting 

2s(t) + s 2 (t) + s 2 (t) = x(t), (56) 

we may write (55) as 

Bf\x(t)\ = my(t), (57) 

where 

f(x) = Vl + x - 1, x 5= -1. (58) 

Then s(t) may be found from the solution x(t) of (57). (This requires 
factoring 1 + x(t) in the form g(t)g(t), where the bandwidth of g is 
half the bandwidth of x.) Then with y = f(x) we have for the inverse 

x = <t>(y) = 2v + y 2 . (59) 

Setting 

00 

x = I m k x k , (60) 

1 

we have from (29) with ai = 2, a 2 = 1, a k = for k 2* 3, 

*i = 2y (61.1) 

* 2 = By 2 (61.2) 

x 3 = -B{y-Hy 2 ) (61.3) 

COMPANDING PROBLEMS 3029 



x 4 = B[yH(y-Hy 2 )] + l/AB(Hy 2 ) 2 (61.4) 

x 5 = -B\y.H[yH(yHy 2 )]\ - l/4B[yH(Hy 2 ) 2 ] 

- l/2B[(Hy 2 )-H(y-Hy 2 )]. (61.5) 

Replacing H by / - B we have from (30) the alternate forms 

X3 = B(yBy 2 ) - By 3 (61.3a) 

x 4 = B[y.B(yBy 2 )) - B(yBy 3 ) + l/4B{By 2 ) 2 

- S/2B(y 2 -By 2 ) + 5/45/ (61.4a) 
x, = Biy-B[y.B(y-By 2 )]) - B[y-B(y-By 3 )) 

- 3/2fi[y-B(y 2 -5y 2 )] + l/4fi[y-(By 2 ) 2 ] 

- S/4B[y-(By 2 ) 2 ] + 5/4B(y-By 4 ) 
+ l/2B[By 2 )-B(y-By 2 )] 

- S/2B[y 2 -B(yBy 2 )] - l/2B[(By) 2 -(By 3 )] 
+ 3/2B(y 2 .By 3 ) 

+ 5/2B(y 3 -By 2 ) - l/4By 5 . (61.5a) 

The factoring of 1 + x can be avoided by developing a series solution 
for s, Bs = s. We have 

x = 2s + s 2 + s 2 



x = mx\ + m 2 x 2 + m 3 x 3 + 



Then setting 



we have 



s = msi + m 2 s 2 + m 3 s 3 + • • • (62) 

s = msi + m 2 s 2 + m% + • • • , (63) 

s 2 = m 2 s f + m 32 Sl s 2 + m 4 {2sis 3 + s 2 2 ) 

+ m 5 (2sis 4 + 2S2S3) + • • • (64) 

s 2 = m 2 s\ + m 3 2sis 2 + m 4 (2s\s3 + si) 

+ m 5 (2sis 4 + 2s 2 s 3 ) + • • • . (65) 

Note that if s belongs to ^(X), then the Fourier transform* of (s + 
is) vanishes outside [0, X] and that of its complex conjugate (s - is) 



♦Here the Fourier transform of s is -i(sgnw)S(o)) where S(w) is the Fourier transform 
of s. 
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vanishes outside [-X, 0]. Thus the Fourier transform of e~ ixt/2 (s + is) 
vanishes outside [-(X/2), A/2]. It follows that the Fourier transform 
of s 2 + s 2 vanishes outside [-X, X], and hence that the sums of the 
coefficients of m n in (64) and (65) are functions whose Fourier trans- 
forms vanish outside [-X, X]. 

Now we can solve successively for s k . It is convenient to introduce 
the Hilbert transform (Quadrature) operator Q 

g = Qg (66) 

to indicate the "hat" of complicated expressions. 
Equating coefficients of m* in 

x = 2s + s 2 + s 2 , 

we have 

s, = Xi/2 = v (67.1) 

s 2 = 1/2x2 - l/2(s? + s?) 

= l/2By 2 - l/2(y 2 + y 2 ), (67.2) 

which may be written, using Bs 2 = s 2 , as 

s 2 = - l/2£y 2 (67.2a) 

S 3 = 1/2X 3 - (SiS 2 + SiS 2 ) 

= - l/2B(y-Hy 2 ) + l/2yBy 2 + l/2y-QBy 2 . (67.3) 

Here we may write 

yBy 2 =yy 2 -y.Hy 2 
and then use Bs 3 = s 3 to obtain 

s 3 = - l/2B(y-Hy 2 ) - l/2B(y-Hy 2 ) + l/2B(y-y 2 ) + l/2B(y.QBy 2 ) 
= - l/2B[yH(y 2 + y 2 )] + l/2fi(y-y 2 ) + l/2B(y.QBy 2 ). 
Then since H(y 2 + y 2 ) = 0, we have 

s 3 = 1/2 B(y-y 2 ) + 1/2 B(y-QBy 2 ) (67.3a) 

s 4 = 1/2 x 4 - (8is a + 1/2 s 2 2 ) - (s^ + 1/2 s|) 
= 1/2 B[y-H(y.Hy 2 )] + 1/4 B(Hy 2 ) 2 

- 1/2 yB(y.y 2 ) - 1/2 yQB(y.y 2 ) 

- 1/8 (By 2 ) 2 - 1/8 (QBy 2 ) 2 . (67.4) 

There appears to be no simplification here. One may prefer the 
alternate expression (61.4a) for x 4 to eliminate the H operator. Note 
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that QB can be replaced by the bandlimiting quadrature operator B 
where 

P" , v 1 — cos X(£ — s) , ,_ ov 

£* = (Bg)(t) = J _ g(s) ^^ ds. (68) 

VII. THE COMPLEX-VALUED COMPANDING PROBLEM 

The c.v. companding problem is considerably more complicated than 
the r.v. companding problem, even for the same analytic companding 
function. For example, if 

f{x) = tan" 1 * + ex, (« > 0), 

we know from the Landau-Miranker theory that the r.v. companding 
problem 

Bf(x)=y 

has a unique solution x in $l 2 corresponding to every real-valued y in 
^ 2 - However, in the case of complex-valued y, this may not be true 
because x must then take complex values which, if the norm of y is 
not restricted, may be singularities of /. In addition, we are confronted 
with the problem of establishing the uniqueness of the solution, which 
may require still more severe restrictions on the norm of y. 

Beurling's uniqueness proof (see Ref. 1) for the r.v. companding 
problem is elegant and simple: Suppose f(x) = j0(\x\), \x\ -> 0, and 
is monotone increasing, and further that 

Bf(x x ) = y and Bf(x 2 ) = y, 

with Xi, x 2 , (and y) in & 2 . Then f(xi) and f(x 2 ) belong to L 2 and B{f- 
(xi) - f(x 2 )\ = 0, i.e., the Fourier transform of \f(xi) - f(x 2 )\ vanishes 
over (-X, X), and therefore (/(xO - f(x 2 )} must be orthogonal to 
(*i - x 2 ). But this is impossible unless Xi = x 2 , for otherwise (x x - x 2 ) • 
{/(*i) - /(*2)1, which is everywhere non-negative, will be positive 
everywhere on the real axis, except at the isolated zeros (if any) of 

(*i - x 2 ). 

For establishing uniqueness in the c.v. companding problem, it 
would seem that the weakest analogue of monotocity should be 
"schlichtness" of /, i.e., that x should be confined to a region, where 
f( Xl ) = f(x 2 ) implies x x = x 2 . This suffices to establish uniqueness of 
the solution in the special case where x has a one-sided Fourier 
transform, but we are not able to see that it suffices in the general 
case. We can establish the following analogue of Buerling's theorem, 
where, without loss of generality, we assume /'(0) is positive. 
Theorem 1: Suppose /(0) = 0, /'(0) > 0, and f(z) is analytic in a convex 
region G including the origin, wherein 
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Re{f'(z)\>0. 

If Xi(t) and x 2 (t) belong to i^ 2 and are confined to G for all real t, then 

Bf(x 1 )=y and Bf(x 2 )=y 

imply 

xi(t) - x 2 (t). 

Proof of Theorem 1: Since G is convex, any two points Xi and x 2 in G 
can be connected by a straight line segment in G. Suppose {x\ — x 2 ) = 
re ,e , where r > 0. Then integrating /' along the connecting line 
segment, we have 

/(*i) " f(x 2 ) = e ie I f'(x 2 + se ie )ds, 
Jo 

and hence 

Be fM ~ /fe) = Re i P /'(x, + S e»)<fe > 0. 

Xi — x 2 r Jo 

In case x x —*x 2 , the limit is Re f'(x 2 ) > 0. 

Now \f[xi(t)] - f[x 2 (t)]\ belongs to L 2 and must be orthogonal to all 
members of & 2 , in particular to \xi(t) — x 2 (t)\; i.e., setting 



Pit) = (*, - x 2 ){f( Xl ) - f(x 2 )\ = |xi - x 2 \ 



2 l/(«i) - /te)} 

*1 - 3:2 



the integral of P(£) must vanish. However, we see that the real part of 
P(t) is non-negative everywhere on the real axis, and vanishes only 
where X\ = x 2 . Since the integral of P is zero, its real part vanishes 
a.e. Thus the function [xi(t) — x 2 (t)\ in <% 2 vanishes a.e., and hence 
everywhere. □ 

Now we can establish that the c.v. companding problem, 

Bf(x) = y, 

will have a solution x, which will take values in a disk centered on the 
origin, wherein Re(/'(z)| > 0, provided \\y\\ 2 is sufficiently small. Then 
the uniqueness of the solution follows from Theorem 1. 

An objectionable, but inherent, feature of companding problems (as 
formulated here) is that a restriction on || y IU is not sufficient to give 
a corresponding restriction on ||jc||«. We can, however, establish that 
|| x || 2 will be small if ||y|| 2 is small, and hence that ||x||„ will be small, 
according to the inequality (given in the introduction) for a function 
g in ^j(X), 

U«iL<C JkRUU. (69) 
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In the sequel, we assume, for convenience and without loss of 
generality, that/' (0) = 1, 

f(z) = z + l b k z\ \z\<R , (70) 

2 

where R Q (perhaps oo) is the radius of convergence of the series. We 
exclude the trivial case f(z) = z, and define 

M(r) - max | f'{z) - 1 1 « 2 k\ b k \ r k ~\ r < R Q , (71) 

l«-r| 2 

which increases steadily from to oo, allowing us to define p uniquely 
by 

M(p) = 1. (72) 

Then it is clear that 

Re{/'( 2 )}>0 for \z\ < p. (73) 

We are now able to establish the following result. 
Theorem 2: Let y(t) be any complex-valued function in & 2 = ^(X), 
satisfying 

Vx7»||y|| 2 =5 max{r[l - M(r)]\ = r [l - M(r )]. (74) 

0<r<p 

Then the companding problem 

Bf(x)=y 

has a unique solution x = x(t) in & 2 - 

Proof of Theorem 2: We can use the method of Landau and Miranker 
to obtain a Cauchy sequence \x n \ converging to the solution x, provided 
we restrict ||y|| 2 , in the end, to be sufficiently small that all the 
approximants satisfy \x n \ < p. 
Assuming the norm of y to be sufficiently small, we take 

x x =y = By, (75) 

which should be a good approximation to x for small y. Then we set, 
following Landau and Miranker, 

x n+ i = x n + y - Bf(x n ), n>l t (76) 

so that, by induction, Bx n = x n , i.e., x n belongs to ^ 2 . We have, writing 
the same equation for n - 1 and subtracting, 

X„ + i - X n = X n - Bf(x n ) - \Xn-\ ~ Bf(x n -i)} 

= B[x n - f(x n ) - {*„-! - /(*„-!)}]. (77) 
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Now we write 

/(*„) - Xn - {/(*„_i) - X„-l) = I {/'(z) - ljcfe. (78) 

Then assuming that 

| x n | ^ r < p for all n ^ 1 (79) 

[and all £, suppressed in the notation x n = x n (t)] we have in (78) 



f 



{/'(*) - U<fe 



M(r)|x„-o: n - 1 | s (80) 



where 

M(r) < M{p) = 1. 

Substituting in (77) the inequality (80) for the magnitude of the 
function in (78), we obtain 

IUn+i - x n \\ 2 ^ M(r) || x„ - x„_i|| 2 . (81) 

So, under the assumption (79), \x n ] forms a Cauchy sequence converg- 
ing to x in ^ 2 [cf. Landau, Ref. 1]. It follows from (76) that 

Bf(x) = y. (82) 

Now we would like to see how large ||y|| 2 may be in order that (79) 
hold, giving the conclusion in (82). We write 

X n = Xi + (X 2 ~ Xi) + (X 3 ~ X 2 ) + • • • + (X n ~ Xn-i) (83) 

from which follows 

»i 

Iknlla^E \\x k -x k -A 2 , (84) 

i 

where 

Xi = v and x = 0. 

Applying (81) to (84), we have 

Unh'S-^-Wyh, (85) 

1 — a 

where a = M(r) < 1, provided (79) holds. This will be the case, 
according to (69), if 

||x„|| 2 ^ yfirj\r for all n > 1, (86) 

which, in turn, will hold if in (85) we have 
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v£7x||yll2*r[l-M(r)]. (87) 

Here we are free to take the maximum over r. Thus the problem 
will have a solution x satisfying the hypotheses, provided the norm of 
y satisfies 

yf^/X \y\%< max{r(l - M(r)\ = r [l - M(r )], (88) 

(Xr<p 

and the solution is unique according to Theorem 1. □ 

We note that if (88) is satisfied, then the solution x satisfies, 
according to (85) and (69), 

|*|. < fb < p. (89) 

So, in fact, the restriction (88) on the norm of y is too severe. We 
obtain a slightly better result later, using a different method. 

We now wish to show that if y\ and y 2 are close to each other, then 
the corresponding solutions, Xi and x 2 , are also close to each other. 
Lemma: Let y x and y 2 satisfy the hypotheses on y in Theorem 2. Then 
the solutions of 

S/(*i) = y\ and Bf{x 2 ) = y 2 
satisfy 

■■ ii I — 7T yi ~ y2 2 _ n mi\ 

l 1 '^ 11 -*^^^* 2 "' (91) 

Proof of the Lemma: We have 

x l -x 2 =y 1 -y 2 - B[f{x x ) - *i - /(* 2 ) + x 2 ], (92) 

giving 

ll*i -* 2 || 2 *s II yi -y 2 || 2 + II5[-]II 2 ^ II yi -y 2 || 2 + IIHIU- (93) 
Also, since, according to (89), 

|*i|^r <p and |x 2 |sSr <p, 
we have from (78) and (80), 

Wf(xi) - xi - f(x 2 ) + X2II2 < M(r ) ||*i - X2II2, (94) 

which with (93) gives 

. . llyi-y 2 || 2 ™ 

*1 — *2 II 2 * - .,, x • V»OJ 

1 - M(r ) 
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This, with (69) and the assumptions on y x and y 2 , establishes the 
lemma. □ 
With this Lemma and Theorem 2 we can show that the problem 

Bf\x(t;m)\ =my(t), 

for fixed y in $l 2 , has a unique solution in @ 2 for all complex m of 
sufficiently small magnitude, the solution x(t; m) being a continuous 
function of the complex variable m in a certain disk centered on the 
origin. To establish for each t that F(m;t) = x(t;m) is an analytic 
function of m in that disk, we show that F(m) has a derivative 
(nondirectional) there. Working with the derivative we are able to 
improve on Theorem 2. It is convenient now to set VA/*Jy|i = 1 so 
that|y(t)|«l. 

Theorem 3: Let y{t) be any complex-valued function in 3) 2 ■ ^ 2 (X) 
satisfying 

VxAII;yll2 = i. 

Then the problem 

Bf\x(t;m)\ = my(t) 
has a unique solution x(t; m) in & 2 for all complex m satisfying 



| m| *£a(p) = I [1 -M(r)]dr, 



(96) 



where M(r) and p are defined in (71) and (72). Furthermore, for each 
fixed real t, x(t;m) is an analytic function of m, \m\ ^ a(p), and hence, 
since x(t; 0) ■ 0, 

x(t;m) = X m k x k (t), \m\ s= a(p), (-00 < t < oo) (97) 
i 

where the x k (t) depend only on y(t) and f. 

We note, before proving Theorem 3, that in Theorem" 2, < r < p, 
and in Theorem 3 

a(p) =1 [1 - M(r)]dr + I [1 - M(r)]dr, 

JO Jr 

where M(r) increases from to 1 over (0, p). Thus 

J [1 - M(r)]dr > r [l - M(r )], 
•Jo 

and hence 

r [l - M(r )] < a(p) < p. (98) 

COMPANDING PROBLEMS 3037 



Then, according to Theorem 3, the c.v. companding problem 

Bf(x) = y 

has solutions for y of larger norm than specified in Theorem 2. 
Proof of Theorem 3: We first consider the solutions x x = x{t;mi), x 2 = 
x(t;m 2 ), corresponding to y x = m x y(t), y 2 = rr^yit), where 

m 2 = m l + e, e=\e\e w (99) 

and 

I mi | + | el < ro[l - M(r )], (100) 

so that Theorem 2 and the Lemma apply. 
We have 

y 2 -yi = ey(t), (101) 

and hence, from the Lemma, 

■*-*fc*T^'*-£S- (102) 

Now 

B{f(x 2 ) - f( Xl )\ = ey, (103) 

which we rewrite as 

*2-*i m _ B \ fM ~ f(xi) - (x 2 - xS [ (1Q4) 

We intend to let c —> (with any argument) and show that the quantity 
on the left tends to a limit, independent of arg(e); viz., F'imtft), where 

F >( m;t ) = A x(t;m), | m | < r [l - M(r )]. (105) 

dm 

From (78), (80), and (102) we have 

II f(x 2 ) - f{Xi) - (x 2 - x x ) h < M(r )|| x 2 - X! || 2 

* i *j vs sy ■ < io6) 

1 - M(r ) 

So 

/<*) - /(x t ) - fa - *i) belongsto L2 . (107) 

e 

We also have from (102), or the Lemma, 
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ll*2-*ilU« - ' ' . . (108) 

1 - M(r ) 

Thus we may write 

x 2 = x 1 + eu, (109) 

where 

u = u(t;m lt e) belongs to _^ 2 and | u \ 

= ^(1) as «-*0. (110) 
Then 

rt*2)-/(*i)-/(Xi + eu) -f( Xl ) 

= f u/'(x 1 )+^V). 

So (104) may be rewritten as 

u = y- B\uf'( Xl ) -u+ j?Uu% (111) 

Now letting e -> and replacing mi by m and x x by x(t;m), we obtain, 
setting u(t;m, o) = u(t;m), the equation 

u(t;m) = y(t) - B\u(t;m)[f ' {x(t;m)\ - 1]}, (112) 

\m\ **A<r [l -M(r )]. 

Here we make the identification 

u(t;m) = F'(m;t) = — - x(t;m), \ m \ =S A (113) 

dm 

by verifying that (112) has a solution u(t;m) in & 2 , in fact, for \m\ 
larger than r [l - M(r )]. We observe, since x(t;0) = 0, and/'(0) = 1, 
that 

u(t;0) = y{t). (114) 

Actually, we can obtain better estimates for | x(t;m) \ by integrating 
its partial derivative from to m. 
We consider the equation for u, 

u = y-B\u.[f'(x)- 1)\, (115) 

assuming x = x(t; m) is known and satisfies 

||x|U«r<p, (116) 

so that 

\f'(x)- 1| < M(r) < M(p) - 1. (117) 

Using this inequality in (115) we obtain 
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lulls < Ilyll2 + M(r)||u|| 2 , (118) 

or 

II -.11 J«r/\ 

(119) 



" ll2 ^l-M(r) = l-M(r)- 



The last inequality implies that (115) has a solution u in 02 
(obtained iteratively), in fact, for || x ||. < p, since M{r) < 1 for r < p. 
We note further that the inequality (119) is crude, with equality 
possible only for r = 0, for we cannot have 

| f'\x(t;m)\ - 1 1 - M(r), (-» < * < ») 

unless xU;m) ■ 0. Therefore, in (115) we have 

\\B\u[f'(x) - l]||2<M(r)||u|| 2 for 0<|x|.«r<p. 

So we have strict inequality in (119) for < r < p. Hence, 

ll"U-<; — 77TT for 0<|*k«r<A ( 12 °) 

1 — Az(r) 

Now let us set 

m = ae», a > (121) 

and 

r(«) = max ||x(t;ae*)|U (122) 

We want to see how large we can make a, say a(p), and have r(a) < 
p. Using (120) in 

*(£;m)= I u(ttf)d|, (123) 



(124) 

Writ)] 



we obtain the inequality 
Then, after defining 

■w-fHfcgi' °«« *«<*>• (125) 

it is clear, since M(r) is an increasing function of r, that we will have 

r(a) < 8(a), < a *S a(p). (126) 

Differentiating (125) with respect to a, we obtain the simple equa- 
tion 
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s'(a){l-M[s(a)]\ = l, 
or, considering a as a function of s, 

da 

— = 1 - M(s). 
Thus 

a(s) = I [1 - M(r)]dr, < s «£ p. (127) 

Jo 

We have s(a(p)) = p and r(a) < s(a) for a > 0. So we will have 

||xtero)|L<p for (128) 

I m | < a(p) = I [1 - M(r)]dr. (129) 

According to (120) and (128), the partial derivative u(t;m) will exist 
for \m\ somewhat larger than a(p). This completes the proof of 
Theorem 3. □ 

VIII. AN ILLUSTRATIVE EXAMPLE 

It will be shown in a future paper that the r.v. companding problem 

B {rriJ* a,y ' x<1 » x,y in M\), (130) 

is equivalent to finding the reproducing kernel for a certain Hilbert 
space of bandlimited functions. The specific problem with X = 2 (for 
convenience) and 

sin It 

is quite easily solved. For real m > —2, the solution is 

/x x nn sin 2t „ 9 /sin t\ 
x(t;m) = 2(3 -^- - P 2 [—J , (132) 

where 

|8 = m/(2 + m). 

We need not be concerned here with the derivation of this solution, 
as we will later show directly that it satisfies 

_ f x(t;m) j sin 2t 

[1 - x(t;m)) 2t 
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for all m in a certain region of the complex plane, but for no other m. 
We know from Theorem 3 that (133) has a unique solution for 
sufficiently small | m\, and that the solution is an analytic function 
of m. It follows that (132) is the solution of (133) for all complex m, 
| m | < | mo |, for some | mo| > 0. Since m = -2 is the only point where 
x(t;m) is not analytic, we might suppose | mo| =2. The series expan- 
sion of x{t,m) certainly converges (uniformly in t) for all 
| m | < 2, but it is not a solution of (133) for all such m. For example, 
we have 

and 

20 



m = 



1-/3' 



Then for = ± itt/2, we have 



m = -, and x\± ^;m ) - I .. 

1 ± i - 

TV 

Therefore, the meromorphic function of t, 

x{t;m) 



f\x(t;m)\ = 



1 - x(t;m) 



will have poles at t = ± tt/2 for m = -2/(1 ± ;2/tt). Thus we have here 
an example, \m\ < 2, for which (132) is not a solution of (133). 
However, it is, according to Theorem 3, for all m satisfying | m \ ^ 
3 - 2>/2~. This, as it turns out, is an overly conservative upper bound 
for | m | . 

We now turn to the problem of determining precisely those m for 
which (132) is a solution of (133). 

First we can easily show that the r.v. problem has no solution for 
m *S -2 by convolving both sides of (133) with l/irK(t), where K(t) = 
(sin tf/t 2 . The result is 



£ 



x{s;m) 1 . , m 



- K(t -s)ds = - K(t). (134) 



Loo 1 — x(s;m) tv 2 

Since x/(l - x) > -1, and K(t) & 0, we have 

f * ( * m) - i Kit - s)ds > - - f K(t)dt = -1 

J-* 1 - xis;m) it it » / -°° 
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This gives, setting t = in (134), 

m > -2. (135) 

Now to proceed towards our stated goal, we first write 

1 - x(f, m ) = (l - fie" 8 M (l - He" »M. (136) 



Then in order for x(t;m) to be a solution of the problem, the Fourier 
transform of the function 

x( t;m) sin It 

h{t\m) = - — — - m — — 137) 

1 - x(t;m) 2t 

must vanish over (—2,2). With a bit of manipulation we arrive at the 
expression 

h(t;m) = g(t;m) + g(-t;m), (138) 



where 



Ml - e» 8in fc 

g(t;m) = YT~ j ^ ( i 89) 

lit ( 1 



-**^) 



and 



9- (m/2) ^i. 

M 1 + (m/2) 
We now introduce the complex variable t = t + iu, and observe that 



| JW + * I1 0|_ J ^-_.J > u ^ + co. (140) 

Then if the denominator satisfies the condition, 

1 - 0e iT ^LI) is zero-free for u > 0, (141) 



it is easy to see (by contour integration in the upper half -plane) that 
G(co;m) = g(ir,m)e- iwt dt = for w < 2. (142) 

On the other hand, if the function in (141) has zeros r k in the upper 
half plane u > 0, (it must have no real zeros in order for g to have a 
Fourier transform) we will have, by the calculus of residues, 
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G(o);m) = £ c*e- ,w *, c fc ¥> 0, (143) 



i 



where n, depending on ft is finite, since it is clear that the function in 
(141) can have only a finite number of zeros in the upper half -plane. 
Since the Fourier transform of h(t;m) is given by 

H(«, m) = G(u;m) + G(-a>;m) (144) 

it will vanish over (-2,2) if, and only if, the condition (141) is satisfied. 
Now the values taken by e ir sin t/t in the upper half-plane, u 5= 
are precisely those values on the boundary and interior of the cardioid- 
like region whose boundary is described parametrically by 

i, sin t 

e lt , -T «S t =£ 7T. 

t 

(Some values are taken more than once.) Then x(t;m) will be a solution 
to the problem except for those values of m such that 1//3 is a point 
on the boundary or in the interior of the cardioid-like region. By the 
mapping 



m = 



»-■ 



x(t;m) is a solution to the problem for precisely those (finite) m lying 
in the region to the right of the boundary line described parametrically 
by 

m = r^ , -TT^t^w. (145) 

t 

This region (see Fig. 1) includes the half-plane Re{m\ 2= -4/3, its 
boundary indenting more to the left near the real axis, having a cusp 
at its leftmost point, m = -2, where it is tangent to the real axis. It is 
found that the distance from the origin to the boundary is minimal 
(see circle in Fig. 1) at the point m and its conjugate, where 

m = (-2 + £) + & 

£ = - = .4895273114 = (2.42786943)" 1 , (146) 

to is the smallest positive root of sin t/t = cos t + sin t, 

I ttioI = = 1.58781760, arg|m } = 7 + to. 

sin t 4 

So I mo I is the largest number such that x(t;m) is a solution for all m 
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Fig. 1— Open region in m-plane (unshaded) for which eq. (132) is a solution of eq. 
(133). Shown is the largest disk (centered on the origin) contained in the region. 



satisfying \m\ < \m \. Also, 

| ra | = min Ro(t), 

t 

where Ro(t) is the radius of convergence of the series 
x(t;m) sin 2t 



= m 



+ m z h 2 (t) + m a h 3 {t) + • 



1 - x(t;m) 2t 

the minimum occurring for t = ±t . 

IX. CONCLUSION 

The expressions for the nth order components x n {t) of the series 
solution to the companding problem become so complicated that, for 
practical purposes, only the first few are of interest. These should be 
useful in correcting small distortions in nonlinear transmission sys- 
tems which fit the companding model. It would appear that the 
corrections applied internally to the inverse function (the z k in Section 
V) would be more effective for correcting larger distortions, especially 
if the lower frequency components are predominant in the signals. In 
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this connection, the simpler "approximate identity" should be quite 
effective for correcting small to moderate distortions of a more general 
nature, as evidenced by the inequality given in the appendix. Experi- 
mental evidence of the effectiveness of these correction schemes would 
be desirable. 

The question of the convergence of the series solution is a matter 
of little practical concern, but the fact that it does attaches more 
mathematical significance to the results. To settle this question we 
had to show that the complex-valued companding problem has solu- 
tions for functions of sufficiently small norm. This generalizes the 
result for functions of one-sided spectra; and whether or not the 
general result will ever find practical application, it is an interesting 
addition to a theory, though still incomplete in many respects. For 
example, it is doubtful that the condition that x(t) be confined to a 
convex region G where Re{ /'} > is a necessary condition for unique- 
ness of the solution. In this connection, one could probably use analytic 
continuation arguments to show that the specific problem examined 
in Section VIII has solutions only for those values of m for which the 
(particular) solution given is the only solution, this being unique for 
sufficiently small \m\, and being an analytic function of m having no 
branch points. Also there is the difficult question of determining for 
what y(t ) the companding problem has a solution, where particular 
interest is attached to the real-valued problem with analytic compand- 
ing functions. It can be shown for the case fix) = x/(l - x), x < 1, 
that the problem has a solution for every (r.v.) y in ^> satisfying 
y > -l. This suggests (conjecture) that the r.v. companding problem 
with f{x) = x/(l - x 2 ), -Kx<l, has a solution for every (r.v.) y in 
&2, or more generally for monotone f(x) defined on (-1, 1) having 
singularities as strong as poles at ±1. In general, it is not enough for 
fix) to increase from — » to +°° over its range of definition in order to 
draw the same conclusion; e.g., fix) = log(l + x),x> -1. The questions 
raised here are certainly worthy of future consideration. 

In connection with the series solution, one naturally inquires 
whether an explicit formula (albeit complicated and involving parti- 
tions of various kinds) can be given for the general term x n . Perhaps 
the combinatorics experts will consider this question. 

We note that the solution x = B<piy), valid for y (of sufficiently 
small norm) whose Fourier transforms vanish outside [0, X], is verified 
by the fact that in (29.n) the expression reduces to x n = a n By n , the 
other terms vanishing because B is operating on functions whose 
Fourier transforms vanish over (-<», X). The same reduction occurs in 
the expression (34.n), because, in this case, B is operating on functions 
whose Fourier transforms vanish over (-°o, 0) and agree over [0, X]; 
i.e., 

3046 THE BELL SYSTEM TECHNICAL JOURNAL, DECEMBER 1983 



B[...] = By n 

holds for each term in (34. n), the sum of all the coefficients being a n . 
Some interesting identities are obtained by equating the expressions 
for x n in the series solution of the general problem to those obtained 
from the explicit solution (given in the introduction) to the special 
problem 

B log(l + x) = my, 

which involve y, the Hilbert transform of y, which does not appear in 
the more general expressions. For example, we find from (34.3) that 
x 3 in the series solution of this problem is given by 

x 3 =l/2B(y.By 2 )-l/3By\ 

and, from the series expansion of the explicit solution, by 

*a = l/SyB(y 2 - y 2 ) - l/QB(yy 2 ) + l/24By 3 + l/4yB(yy). 

It is an interesting exercise to show directly that these two expressions 
are identical. 

Finally, since truly bandlimited signals exist only as mathematical 
abstractions, some attention should be given to developing a mathe- 
matical theory of practical companding problems, 



£ 



f\x(s)\k(t - s)ds = y(t), 



where k(t) is the (absolutely-integrable) impulse response of a practical 
low-pass filter, so that the theory may be extended to signals that are 
merely bounded. Here one may not be interested, for various reasons, 
in the exact solution of this problem but, instead, a compromise 
problem, where the equation is nearly satisfied with both x(t) and y(t) 
being close to bandlimited functions. For example, in many cases 
f\x(t)\ is given (say) by an nth order differential operator acting on 
y(t). Then the (exact) solution x{t) = <p\f(x)\ may be far from a 
bandlimited function. However, if y(t) is close to a bandlimited func- 
tion there should be an approximate solution which is also close to a 
bandlimited function. A case in point is found in Landau's simulation 
of the iterative solution of the companding problem (Ref. 2), where, 
in fact, the equation he was obtaining approximate solutions to was 
the case MO = k(t), (approximately bandlimited) for which the unique 
solution is (a multiple of) the Dirac delta function. 
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APPENDIX 

Suppose f(x) is a monotone increasing function of the real variable 
x, satisfying 

(i) /(0) = 

(ii) < mi < f'(x) *S m 2 < °°, (-« < x < »). 

Then / has an inverse 

(Hi) x = 4>[f(x)\, (-oo<x<oo), 

satisfying, since 1 = f'(x)<t>'\f(x)\, 

(iv) 0< — < <b'(y) ^ — < oo, (-oo<y<oo). 

m 2 mi 

Now let x = x(t) be any function in ^>. We wish to establish 

U-B<f,\Bf(x)}U*y\\xh (147) 

where 7 = An T •»> e ~ ~ ~ x - 

4(1 + €j mi 



Set 
Then 
Now set 



y(t) =y = Bf(x). (148) 

f(x) =y + h, Bh = 0. (149) 



*i = B0(y) (150) 
so that, since x = <j>(y + h) = Bx = B\<j)(y + h)\, 

x-x,= B\4>{y + h) - 0(y)). (151) 
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Since | <j>(y + h) - 0(y)| s= — h, we see that 

mi 



and hence that 



<t>(y + h) - <f>(y)\\ 2 *z — \\h\\ 2 

TTli 



x- x 1 \\ 2 ^— \\h\\ 2 , 

TYl\ 



but we can improve this inequality by writing 

fy+h 

0(y + h) - <f>(y) = I |0'(£) - a \dZ + ah, (152) 

%fy 

where a = - I 1 ). 

2 \rai m 2 



Then, since 



'^'""'^t-i")' < 153 » 

I \rrii m 2 J 



we have 

<t>(y + h) - 0(y) = u + ah, (154) 

where u = u(t) and 

|ul *ife-^) |A| - (155 > 

Thus 

x - xi = B(u + a/i) = Bu, (156) 



and hence 



*-*i|| 2 « |u|a<£(- -) \\h\\ 2 . (157) 

z \mi m2> 



Now we need an inequality of the form || h \\ 2 =S c || x \\ 2 . We have 

h = Hf(x), (158) 

where H = I — B is the high-pass operator. So, clearly 

\\h\\ 2 ^ ||/(jt)|| 2 *Sm 2 ||x||2. 
We can improve this inequality by setting 

v(t) = v= f{x) - ftx, (159) 

where ft = l/2(mi -I- m 2 ). 
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Then 

\v\ * l/2(m 2 - mi)|x|, (160) 

and, since Hx = 0, 

h = Hf(x) = H[f(x) - (3x\ = Hv. (161) 

Hence 

II Jill 2 ^ Bulla *S l/2(m 2 - m 1 )\\x\\ 2 . (162) 

This, with (157) gives 

II x - x x || 2 ^ t( ) (m 2 - mi)||*|| 2 = 7 \\x\\ 2 , (163) 

4 \mi m 2 / 

which is the result (147). The number 7 in the inequality is invariant 
under the interchange of and / in the approximate identity (147), as 
we would expect from using only (ii) and (iv). 

So xi = B4>\Bf(x)\ is a good approximation to x if 7 is small. The 
manipulations leading to the inequality (163) suggest an iterative 
scheme for solving, given y in ^ 2 , 

BF(x) = y, x in tf* (164) 

provided 7 < 1, which will be the case if (m 2 /mi) < 3 + 2 V2. 



We set 




x n = B<f>{y + h n -i), n > 1, 


(165) 


where 




h„ = Hf(xJ, n 3* 0, 


(166) 


and Xo = ho = 0, 




giving x x = B0(y) as in (150). 




Now we wish to show that 




II x - x n \\ 2 «£ 7"ll*ll2, »>1. 


(167) 


We have 




x - x n m B\4>{y + h)- 4>{y + &»-i)}- 


(168) 


Following the previous pattern we write 




<t>(y + h) - <t>(y + h n -i) = u n + a(h - hn-i), 


(169) 


where 




ry+h 

u n = I (*'(*) -«}<& 

*'v+h._i 
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and hence 

\Un\ ^l(—-—) \h - /i„-.,|. (1/0) 

Z \m.i m.2/ 

Then 

x - x n = B\u n + a(h - hn-t)\ = Bu n , (171) 

giving 

II x - x n \\ 2 ** I (— - — ) || h - h n - X \\2. (172) 

# l \m,\ mi) 

Now 

h - hn-i = H{f(x) - fixn-t)}, n&l, with x = ho = 0. (173) 

Here we write 

(*) - f(Xn-l) = J {/'(*) - 0}d$ + P(x - X n -,) 



= V n + P(X - X„-i), (174) 



where 



Tfio — TtX\ 
\V n \ ^ " \X-Xn-A- (175) 

Then 

h - /i n _! = H\u n + 0(x - ac-i)) = Hu n , (176) 

giving, with (175), 

11^ - fc»-ill2 *£ 2 II x - Xts-ills- (177) 

This, with (172), gives 

II* - Xn\\ 2 < tII x - Xn-ilU, (178) 

whence follows, with *o = 0, 

IIx-xJIzsSt'IjL, »>1. (179) 

Note that there is a bonus attached to *i = B0(y), in that only one 
filtering operation is required to obtain it. Thereafter, two filtering 
operations are required to obtain x n from y and x n -\. 
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